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1. Introduction 

QCD at finite chemical potential /i has a complex fermion determinant, which prevents direct 
application of the standard simulation methods of lattice QCD. We adopt the phase-quenched ap- 
proximation, where we replace the fermion determinant by its magnitude, which allows the use of 
standard simulation methods. Outside of the superfluid phase in the {T,}x) plane, it is likely that 
the phase structures of full and phase-quenched QCD are the same. We use exact RHMC simula- 
tions [|l|], employing a speculative lower bound to the spectrum of the Dirac operator to enable the 
required use of rational approximations [^. 

We simulate 3-flavour lattice QCD on 8^ x 4, 12^ x 4 and 16^ x 4 lattices for masses near to 
the critical mass at = 0, for temperatures close to the transition temperature from hadronic 
matter to a quark-gluon plasma. If the critical mass increases with increasing pL, then for masses 
just above nif (0) we would find a critical endpoint at m = mc{li) for small ji. However, as our 
simulations indicate, decreases with increasing pL and no critical endpoint is found for small /i 
1^]. We introduce /i/ = 2/i, which has the interpretation of an isospin chemical potential (at least 
for even numbers of flavours). 

We are now simulating this same 3-fiavour lattice QCD for a range of }X and j8 values outside 
the superfluid phase, in order to calculate the equation-of-state for phase-quenched QCD. 

2. Simulations and results 

The nature of the transition is best determined on finite lattices using the 4th-order Binder 
cumulant for the magnetic order parameter. Since, at finite quark mass, we do not know this order 
parameter, we use the chiral condensate, introducing (hopefully small) finite size effects. The 
Binder cumulant for any observable X is defined by [Q| 



If there is a critical endpoint at small /i/ for m > (0), then the Binder cumulant should decrease 
from its crossover value B4 = 3, passing through the Ising value B4 = 1.604(1) at the endpoint, and 
falling towards its first-order value B4 = 1 as Hj is increased. Small Hj means jJ.] small enough to 
lie outside of the superfluid phase. For T = 0, this means Hj < niji. 

Since, as determined in these simulations, mc{lJLi = 0) = 0.0265(3), we perform simulations 
at quark mass m = 0.02, m = 0.025, m = 0.03 and m = 0.035. At the 3 larger masses we simulate 
at jXi = 0, /I/ = 0.2, jU/ = 0.3, while for m = 0.02 we only simulate at jU/ = 0. On the 12^ x 4 lattice 
where we have the highest statistics, we generate 300,000 trajectories at each of 4 j8 values close 
to j8c (the transition value), and use Ferrenberg-Swendsen reweighting ^ to continue to j3c, taken 
as the j8 which minimizes B4{\j/Y)- Five noisy estimators of per trajectory are used to obtain 
an unbiased estimator for B4. 

Figure |l| shows these Binder cumulants as functions of /X/ for m = 0.03 on 8^ x 4, 12^ x 4 
and 16^ x 4 lattices. Rather than decrease with increasing /i/, the curves for the larger lattices 
show a modest increase with increasing Hj, and hence no sign of a critical endpoint. The graphs for 
m = 0.035 are similar, except that we do not have 'data' on 16^ x 4 lattices. At m = 0.025, although 
there is a suggestion of such an increase, but the 'data' is consistent with no Hi dependence. 
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Figure 1: Binder cumulants at T = Tc as a functions of /x/ at m = 0.030. 
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Figure 2: Binder cumulants ?XT — Tc as a functions of max jii — 0.2. 
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It is instructive to examine the behaviour of S4 as functions of m at fixed /i/ for evidence of 
a critical point in the universality class of the 3-dimensional Ising model. We find that for each of 
our 3 pLi values the curves for the S'' x 4 and 12^ x 4 lattices cross very close to the Ising value. 
Not only does this indicate that the critical point belongs to the Ising universality class, but also 
suggests that using i/Zi/A as the magnetic order parameter is a reasonable choice. Figure § shows 
S4 as functions of m for jXi = 0.2. Estimating the position of the critical points as the values of m 
where the 12^ x 4 curves achieve their Ising value gives: mc(0) = 0.0265(3), nic(0.2) = 0.0259(5) 
and mc(0.3) = 0.0256(4), i.e. mt (jU/) is a slowly decreasing function of /i/. 

From the same 'data' we have also calculated the chiral susceptibility 



(2.2) 



where the i/ai//s on the right-hand side are lattice averaged quantities. Finite size scaling at the 
critical point predicts that 

Xi|,^i,{^T,)=Lix■ (2.3) 

y 

Hence if we plot L^v X(j,xf,{L,Tc) as functions of m for different Ls, the curves should intersect at 
the critical point. In figure ^ we plot this quantity for /X/ = 0.2, and note that curves for 8^ x 4 and 
12^ X 4 lattices would intersect between m = 0.025 and m = 0.03 which is where we found the 
critical point from the Binder cumulants. Similar results obtain for /i/ = and Hj = 0.3. 
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Figure 3: Finite size scaling for the peak of the chiral susceptibilities at ^/ = 0.2 



3. Equation of state 

The equation-of-state (EOS) expresses the pressure p, the entropy density s and the energy 
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density e as functions of temperature T and /X/. (Calculations of the equation of state for QCD at 
finite T and /i have been performed for example in ||^ and [0].) The pressure p is simply related 
to the partition function through: 

p='^\nZ{T,pLi). (3.1) 

However, we do not actually measure the partition function Z in our simulations, only observables. 
Z(r,0) can be calculated by numerically integrating 



d\nZ , V 



(3.2) 



d^ \ J 

where Sg is the plaquette action. We can then numerically integrate 

d\nZ _ NfV 3 

at constant jS, where '^jl is the the isospin density, to obtain Z{T,Hj). Figure ^ shows the Hi 
dependence of at fixed j3 values. For the upper 2 j3s we also have 'data' for larger /i/S, up to 
saturation which occurs for /i/ ^ 2. 



(3.3) 
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Figure 4: Isospin density as functions of jj.] at fixed /3 values 



To obtain T in physical units requires knowledge of the running of the coupling constant, 
(5 = j5{a). This is determined at /i/ = 0. Once this running of the coupling constant is known, this 
can be used to determine £, since 

d 

e = — ^InZ. (3.4) 
V dT ^ ^ 
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4. Discussion and conclusions 



For li small enough or T large enough that phase-quenched QCD is in its normal rather than 
its superfluid phase, full QCD and phase-quenched QCD are expected to have the same phase 
structure. Lattice simulations indicate that the fluctuations of the phase of the fermion determinant, 
on lattices large enough to observe this phase structure, are sufficiently small for this to be so |||]. 
Random matrix calculations in the epsilon regime agree with this conclusion 

For 3-flavour lattice QCD our phase-quenched simulations show no evidence for a critical 
endpoint in the range of m and pL values where it would have been expected if it were associated 
with the critical point observed when the quark mass is varied. This is in agreement with the results 
of de Forcrand and Philipsen obtained using analytic continuation methods from simulations at 
imaginary jx, for full QCD [|T0|]. 

We suggest that the softening of the transition as jXj is increased is because the introduction 
of an isospin chemical potential reduces the symmetry (at least for even numbers of flavours). 
Reducing symmetry tends to soften transitions. For example, reducing the number of flavours from 
3 to 2, reduces the symmetry at m = 0, = from SU{'i) x SU{1>) to SU{1) x SU {!) and the 
transition softens from first order to second. The addition of a mass breaks chiral symmetry to a 
pure vector symmetry, and softens phase transitions to crossovers. 

Determination of the true scaling field for the magnetic order parameter - a linear combination 
of the chiral condensate, the plaquette action and the isospin density - is needed to remove the 
finite size effects. Use of the methods of de Forcrand and Philipsen to improve the signal/noise in 
calculating the jx} dependence of the Binder cumulants, would also help [11, 12]. 

Fodor and Katz have predicted a critical endpoint at Te = 162(2) MeV and /i = 120(13) MeV, 
using their reweighting methods [13]. However, this /i is considerably beyond mji/2 and thus 
beyond the reach of both phase-quenched and analytic-continuation methods. 

We are performing simulations along lines of constant j8 to determine the equation-of-state for 
phase-quenched QCD to compare with that for full QCD [j^ ^ [T^]. At low jSs we are restricted to 
/I/ < rrij,;. At high jSs, where the system is in the plasma phase for all ^U/S, we can cover the whole 
range of jXi. What can we learn from the resonance gas model [|l^] and chiral perturbation theory 
about this equation-of-state? 
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